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Transistor Surface Effects Responsible for 

Anomalous Third-Order Intermodulation 
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Telephone System 
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(Manuscript received June 12, 1979) 

Anomalous intermodulation distortion has been observed in the 
SG submarine cable repeater amplifier. The anomalous distortion 
was found to be dependent upon the crystal orientation in the SG 
bipolar transistor. This fact suggested that transistor surface effects 
might be responsible for the anomalous distortion. To test this hy- 
pothesis, computer models of the transistor with surface effects were 
used to compute the third-order modulation coefficient Ma. Surface 
effects in the collector- base junction alone were modeled by a junc- 
tion-controlled mos capacitor with fast surface states. The mos ca- 
pacitor model then was added to a modified Gummel-Poon transistor 
model to simulate surface effects in the complete transistor. It was 
found that intermodulation contributed by the surface effects was 
dominant in the case of (111) orientation, and the behavior of the 
computed intermodulation was similar to the observed behavior. This 
paper describes the transistor modeling and the intermodulation 
computations developed in testing the hypothesis. The principal con- 
clusion confirms that the anomalous distortion from the SG transistor 
is due to surface effects, in particular, to a high surface-state density. 

I. INTRODUCTION 

The deviation from perfect linearity in a repeater amplifier for a 
high-capacity, long-haul system is very small by normal standards. 
However, when the signal passes through many repeaters, the effect of 
nonlinear distortion can accumulate to a significant value. The result- 
ing modulation noise of many signals is similar to white noise and is an 
important factor in the specification of such system parameters as 
repeater spacing and power-handling capability. Permissible third- 
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order harmonic distortion at the output of a repeater, for a milliwatt 
of fundamental, is typically 10" 10 to 10" 13 mW. 

Measurements of the third-order intermodulation product for the 
SG submarine cable amplifier showed unusual behavior as power was 
increased. 1 The deviation from the expected 3-dB increase in third- 
order power for a 1-dB increase in the three fundamental powers was 
largest for the case where the fundamental frequencies were slightly 
lower than the product frequency, as illustrated in Fig. 1. Distortion 
measurements on the open loop amplifier produced similar results, so 
the feedback network was not a cause. When the bias voltage and 
current were varied separately for each of the signal-carrying bipolar 
transistors in the amplifier, the collector voltage nonlinearity showed 
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Pig. 1— Distortion in the closed loop SG amplifier. 
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the most unusual behavior. Two types of silicon transistors with 
identical geometry but different epitaxial crystal orientations, (111) 
and (100), were tested in the amplifier. The transistors with (111) 
orientation showed the unusual distortion characteristics described 
above, while the transistors with (100) orientation showed classical 
performance. The distortion performance of the collector-base junction 
alone showed unusual performance for the < 111) orientation and 
normal performance for the < 100 > orientation, as shown in Fig. 2. 

These measurements suggested that the collector-base junction was 
the prime source of the unusual distortion. The strong dependence of 
modulation coefficient Mz on crystal orientation suggested surface 
effects as the basic cause of the unusual intermodulation, since surface 
charge is known to be dependent upon crystal orientation. The SG 
transistor chip is shown in Fig. 3. The base contact metalization 
overlaps the collector region with an oxide layer for isolation. This 
forms a pn junction-controlled mos capacitor. It is known that the mos 
capacitor across the base-collector junction can produce a degradation 
of transistor distortion. 2 However, the measured junction capacitance 
for the SG transistors was well-behaved and showed none of the 
anomalous ripples as a function of voltage, as is expected from large 
concentrations of surface states. The Mz measurement of the collector- 
base junction is a more sensitive test than capacitance measurement 
and clearly showed anomalous behavior for < 111) junctions. 

This paper shows that the unusual distortion measured for the 
collector-base junction of the SG transistor is accurately modeled for 
both the (111) and (100) crystal orientations by a junction-controlled 
mos capacitor with fast surface states. The mos capacitor model is 
added to a modified Gummel-Poon transistor model, including the 
effects of base pushout, lateral spreading, and internal heating to 
simulate the measured distortion for an SG transistor. 

The principal conclusion of this work has been to confirm the 
hypothesis that the anomalous distortion from the SG submarine cable 
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Fig. 2 — Mz vs reverse bias for the collector-base junctions of SG transistors. 
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Fig. 3 — Transistor chip. 



transistor is due to surface effects, in particular to a high surface-state 
density. Consequently, it is important in designing transistors for 
applications where intermodulation distortion is critical to take special 
pains to control and minimize surface effects. This may involve using 
silicon wafers that are <100> oriented as well as minimizing the 
collector-base overlap capacitance. Or it may involve special processing 
steps such as: 

(i) Routine characterization of transistor surfaces for surface-state 
and fixed charge densities. 

(ii) Special annealing or gettering steps to reduce the number of 
surface states. 

Appendix E contains a list of symbols used throughout the paper. 

II. DESCRIPTION OF COLLECTOR-BASE JUNCTION MODEL 

The transistor is made up of two base regions separated by a base 
pad, as shown in Fig. 3. Insulated by an SiCVSiN layer, the base pad 
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extends over the collector up to the collector-base junction on each 
side, as can be seen in Fig. 4. This forms a junction-controlled mos 
capacitor with the voltage applied to the junction appearing also across 
the mos capacitor, as shown in Fig. 5. 3 

2.1 Nonlinear small-signal ac model 

A rather complicated ac circuit model for a junction-controlled mos 
capacitor 4 is shown in Fig. 6. Cj represents the junction capacitance 
between the p-base and n-collector, C ox the oxide capacitance, and d 
the capacitance between the semiconductor-oxide interface and the n- 
type bulk. C s models charge storage in the fast interface states at the 
oxide-semiconductor interface. The resistances characterize the follow- 
ing physical processes: 

R nB and R ps — The transition of electrons from the conduction band 
and holes from the valence band into surface states. 

RnD and R p d — The flow of electrons and holes through the surface 
space-charge layer. 

R p b — The generation (or recombination) of holes and their flow rate 
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Fig. 4 — Cross section of SG transistor. 
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Fig. 5 — Schematic of junction-controlled mos capacitor. 
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in the region of quasi-neutrality in the bulk adjacent to the space- 
charge layer. 

Rg—The flow of electrons and holes through the bulk of the semi- 
conductor. For an n-type collector, Ri models the average drop in 
voltage from the base to the semiconductor-oxide interface and helps 
define the surface quasi-fermi level for holes \pF P . The inversion charge 
is stored in Ci and minority carriers can flow to the surface from the 
bulk through R p or through Rj. Holes may be trapped by surface states 
through R PB and stored in C m . However, for an n-type collector with 
positive fixed charge in the oxide, the action of the pn junction on the 
mos capacitor is such as to inhibit the semiconductor surface from 
becoming inverted. 

Since electrons are the important carriers for operation in the 
depletion region, the equivalent circuit shown in Fig. 6 has been 
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Fig. 6— Alternating-current circuit model for a junction-controlled mos capacitor on 
n-type semiconductor. 




Fig. 7— Simplified ac circuit model of a junction-controlled mos capacitor. 
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simplified to that in Fig. 7. Electrons are supplied to the surface from 
the bulk through Rd which is given by Ref. 4: 

1 [ d dx 

QflnA J MX) 

where 

/i„ = electron mobility 
no = dc electron concentration 
d = depletion layer width 
A = surface area. 

The R a , C„ network in Fig. 7 is a simple small-signal model for capture 
and emission of electrons by surface states all located at one energy 
level within the bandgap. In reality, they are spread throughout the 
bandgap, and the admittance for a continuum of states of density 
N 8 (E ) in the bandgap is given by Ref. 5: 

. q ( Ec N a (E)f (l-fo) JV 
Y 88 =jw—\ — dE, (2) 

kT h v 1+ ^t 

CnflsO 

where 

/o = [1 + g 8 exp(u — uf)]' 1 is the fermi function at potential u 
for a quasi-fermi potential Uf expressed in units of the 
thermal voltage kT/q 

g s = ground state degeneracy (= 2 for donors, = 1/4 for accep- 
tors) 

c n = electron capture probability 

/iso = electron density at the silicon surface established by the 
bias 
Ev, Ec = energy levels at the edges of the valence and conduction 
bands. 

Some authors have argued that the integrand of (2) peaks sharply 
about the quasi-fermi level with a width of about kT. e This would 
make it easy to integrate (2) if N S (E) and c n do not vary much with E 
over a range of kT. However as is shown in Appendix A, this is the 
case only at low frequency and high semiconductor doping levels to 
which these intermodulation computations are not restricted. To pre- 
serve accuracy, therefore, integrals of the type exemplified by (2) have 
been evaluated numerically in this work. 

In depletion, the time-constant dispersion measured for mos capac- 
itors is much broader than that predicted from (2). Good agreement is 
obtained if one modifies (2) by assuming a statistical fluctuation of 
surface potential in the plane of the interface. 6 As a first attempt to 
study intermodulation due to surface states, the surface admittance is 
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rederived to retain nonlinear terms up to third order for a continuum 
of surface states. The statistical fluctuation of surface states is taken 
into account via correction terms derived in Appendix B. 

Ml. COMPUTER ALGORITHMS FOR INTERMODULATION 
CALCULATIONS 

3.1 Volterra series representation 

We are interested in computing the third-order modulation coeffi- 
cient for the circuit of Fig. 7. A convenient way of doing this is by 
expressing the current through the nonlinear element as a Volterra 
series of the voltage across the element. 7 If the element is memoryless 
(purely resistive), the Volterra series simplifies to a power series: 

i(t) = hiv(t) + h 2 [v(t)f + h 3 [v(t)f + • • • . (3) 

If the element has memory (capacitive), the first three terms of the 
Volterra series are: 



= hi(t-T)v( 
Jo 

-IT 

Jo Jo 
h(t) - I I I hs(t - ti, t - T 2> t - t 3 ) J] v(r t ) dn. 



ii(t)= fait - t)v(t) dr (4) 

Jo 

W) = h 2 (t -r u t- t 2 ) f[ v(n) dn (5) 



(6) 



The h*a in the integrands are known as Volterra kernels of the first, 
second, and third degree. In the Fourier transform domain, these 
relations become 

/ 1 (w)=H 1 ((o)V(w) (7) 

2 

h(ui, "2) = H 2 ((Ou W2) II VM W 

3 

hfai, w 2 , w 3 ) = Hs(ui, u 2 , 0)3) II VM- (9) 

i-l 

The transformed Volterra kernels Hi(u), if 2 (wi, W2), and H 3(0)1, W2, w 3 ) 
may be found from a knowledge of the physics of the device as 
described in the following sections. 

3.2 Volterra kernels for the surface-state admittance In equilibrium 
(nonlinear case) 

The following derivation is strictly valid only for equilibrium con- 
ditions, which is a good approximation for an mos capacitor in deple- 
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tion with no surface inversion layer. The incremental current flowing 
from the conduction band for an n-type semiconductor into donor or 
acceptor surface states within the energy band dE is given by 

di 8 (t) = -{qN s (E)c n [l - f(t)]n s (t) - qN s (E)e n f(t)} dE, (10) 

where 

i s (t) = time-varying current flowing into the surface states 
q = electron charge 
c„ = electron capture probability 
f(t) = fermi distribution function 
n s (t) = surface concentration of electrons 

e„ = electron emission probability from the surface state. 

The incremental current is also equal to the change in charge per unit 
time 

di s (t) = -qN s (E)^fdE. (11) 

at 

Substitute (11) into (10) to get 

^■ = Cn(l-f)n s -enf. (12) 

at 

This is a nonlinear equation for / in terms of n s with memory due to 
the time derivative. Let n„ be the forcing function. Then one desires 
the small signal components of /"about a dc operating point fo for small 
variations in n s about n^. So let 

f=fo + 8f, n s = n !s0 + 8ns. (13) 

Substituting (13) into (12) produces a dc equation and an ac equa- 
tion. 

c„(l - fo)riso .. . .... 

e„ = (dc). (14) 

/o 

^-8f=c n (l-f )8n >i -^8f-c n 8fSn ll (ac). (15) 

at fo 

Since the nonlinearity (third term on the right) has memory, a Volterra 
series may be used to express 8f. 

8f= Hi(u) a ) o 8n s + Hz{oi a , u)b) o 8n'i + Hz(u a , to&, ic c ) o 8nl, (16) 

where the driving signal 8n s is assumed to consist of three terms: 

8n s = a\ cos(o>i£ + <f>i) + a 2 cos(a>2< + <f>2) + 03 cos(a>3* + fa) (17) 

and the circle operator o indicates that the magnitude and phase of 
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each term in 8n 8 are changed by the magnitude and phase of 

H„(tOi, • • • , "«)■ 

The Volterra kernels Hi, H 2 , H 3 are derived in Appendix C. The 
current into the surface states is computed by converting (11) into the 
frequency domain by using the theorem in Appendix D for differen- 
tiating Volterra series and substituting (16) for 5/ to get 

di s (u>) = -qN s (E)[juaHi(wa) o 8n 8 +j(a a + o3 b )H 2 (o3 a , o3 b ) o 8n 2 8 

+ j(o3 a + 0)b + 03 c )H 3 (o3 a , 10b, u c ) o 8n 3 8 ] dE. (18) 

The total current i 8 (u) is found by integrating di s over the band gap, 

or 

is(a) = -ja a q N 8 (E)Hi(o3 a ) dE o 8n 8 
Je v 

J-E c 
N 8 (E)H 2 (o3 a , o3 b ) dE o 8n 2 8 
E v 

J-E c 
N 8 (E)H 3 (o3 a , wb, ojc) dE o 8nl 
E v 

= Gi(03 a ) O 8tls + G 2 {o3 a , 03 b ) O 8nl + G 3 (o3 a , 03 b , 03c) O fi»f . (19) 

In order to use eq. (19), the surface electron concentration 8n s must 
be known. It may be computed as a nonlinear function of the surface 
potential ^ s and the surface quasi-fermi level for electrons \$>Fn to be 

n s = m expty. - «M/Vr]. < 20 > 

Using a Taylor series expansion 

8n 8 = bi(fy. - S^Fn) + b 2 (84> 8 - 6#„) 2 + 6 3 (% - «W 3 , (21) 

where 

. riso bi b 2 

Note that the nonlinearity does not have memory, so a Taylor series 
expansion is sufficient to describe 8n 8 . Since i, = G(n 8 ) and n„ = B(\p 8 
- \pFn), where B and G are functions from (19) and (20), then i 8 = 
G[B(\}/ S - \pFn)] = D(ip 8 - fan), where the composite function D has a 
Volterra series as shown below. 7 

h = DlM O (8xfr 8 - Sfan) + D 2 (03a, 03 b ) O (% - tyjfc) 2 

+ D 3 (o3 a , 03 b , 03 C ) O (ty. ~ 8^ F n)\ (22) 
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where 

Di(io a ) = G^ a )bi (23) 

D 2 (w a , oj 6 ) = Gi(u a + o>b)b 2 + G 2 (o3 a , 03b)b\ (24) 

Da(u a , o3 b , w c ) = Gi(w a + wfc + w c )& 3 



+ 2G2(w a , <o&, + coc) 6162 + G 3 (w a , 03 b , o3 c )b 3 i (25) 
and 

G^Wa, 03b + (Oc) = - [(^(Wa, Wft + 03 c ) + G 2 (o3b, 03 a + 03 c ) 
+ G 2 (03 C , 03 a + 03b)]. 

3.3 Volterra kernels lor the collector-base junction capacitance and 
surface depletion capacitance 

These nonlinear capacitances can be represented as Taylor series 
expansions of the voltage 8Vc across the capacitance. 

C(8Vc) = Co + Ci dVc + CAVl. (26) 

The current i c through the capacitor can be expressed as 

U = j ( (coSVc + \ CxSVl + i C 2 8Vh\ . (27) 

Converting i c to the frequency domain gives the Volterra series de- 
scription 



U e = jJoSVc + i &8V c + i C 2 5F 3 c) . 



(28) 



The depletion capacitance Cd can be found by differentiating the 
expression for charge in the depletion region (34) with respect to the 
surface potential. This procedure leads to the following expressions for 
the depletion capacitance Volterra kernels: 



Hx («) = jo3 ^^- = - juKH (29) 

civ 



where 



H = exp (u B - UFn) - exp (ufp - u a ) + 2 sinh UFn (29a) 

K = 1 Qi7r (29b) 

G = exp (u s - u Fn ) - exp (-u Fn ) + exp (u Fp - u 8 ) - 

exp UFn + 2u a sinh UFn (29c) 
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and 



L D = 


kT € S ' 

q 2qtii 


1/2 


H 2 (wi, W2> = J 


(tOi + W2) 


K 



— — exp(« s - UFn) — exp UFp 
2Cr 



K [3H 3 3H 



(30) 



[exp («Fn - Us) + exp (« 8 - u Fp )] (31) 

+ exp (u s - u F n) - exp (uf p - u a ) \ , 

where Vr is the thermal voltage kT/q. 

3.4 Finding the dc solution for the MOS capacitor 

Finding the dc solution for the mos capacitor is simply a matter of 
finding the potential at the oxide-semiconductor interface correspond- 
ing to a prescribed dc bias voltage. The charges existing in the 
neighborhood of the interface can be expressed in terms of this surface 
(interface) potential. Since the entire mos capacitor is electrically 
neutral, the summation of these charges must vanish. That is, 

Q e + Q s + Q ss + Qo = Ftys) - 0. (32) 

The charge on the gate Q B is given by 

Q e = C ox {V-^Ms-^ a ), (33) 

where V is the applied voltage, Vms the metal-semiconductor work 
function, \p s the surface potential, and C ox the oxide capacitance. 

The charge in the semiconductor (space-charge region) is given by 
Ref. 3: 

Q a = 2qniL D [exp(u s - uf) - exp(-Mf) + exp(a^ —u,-v) 

- exp(u F -v) + 2u s sinh u F ] 1/2 , (34) 

where u 8 , uf, and v are the surface potential, electron quasi-fermi 
potential, and applied voltage expressed in units of the thermal voltage 
kT/q. 

The amount of charge in the surface states depends on the type of 
state (donor or acceptor), the surface state distribution, and the 
occupancy of the state as described by the fermi function. Figure 8 
shows the charge on the surface states. We have made the conventional 
assumption that the states above mid-gap are acceptors and those 
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Fig. 8 — Charge in surface states. 

below mid-gap are donors. 8 ' 9 Acceptors are negative and donors are 
neutral when filled by an electron. For simplicity, we have assumed 
that all states below the quasi-fermi level are occupied by electrons 
and all above are vacant. Thus the charge in the surface states Q ss has 
been approximated by integrating the surface-state density N s between 
uf and u s assuming all these states to be charged. 



(35) 



Qss = q\ Ns du, 

J i 

Q is the fixed charge in the oxide. 

The value of u s which satisfies eq. (32) is the dc solution. 



3.5 Computing circuit intermodulation 

The modulation coefficient M a measured in decibels is defined by 
Ref. 10: 



Ms = 10 log 



Po(tOi + 052 — W3) 

Po(a>i)Po(u 2 )Po(a> 3 ) 



- 15.6, 



(36) 



where Po is the output power measured in milliwatts for an input 
signal of the form 



Ie = I\ COS COit + I2 COS U)2t + h cos lOst. 



(37) 



The perturbation method is used to compute the circuit response at 
the various combined frequencies as outlined below. 

(i) Solve for the first-order response at wj by replacing the nonlinear 
elements in the equivalent circuit (Fig. 7) with the first-order terms of 
the respective Volterra series, that is, C/ o S\^ s for 8i<j, C>o o 8V for 8ij, 
and D\ o 5^ s for 8i s . The small signal response 8V is then computed 
for the resulting linear circuit for the input /icos uit. 

(ii) Solve for the first-order response at W2 and 103 in a similar 
fashion. Note that D\ is frequency-dependent, so the first-order equiv- 
alent circuit is different for each frequency. 
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(Hi) Solve for the second-order response at o>i + W2 by removing the 
I g generator and inserting intermodulation current generators evalu- 
ated from the second-order terms in parallel with the linear-controlled 
current sources due to the first-order terms (see Fig. 9a). For example, 
for Cj an intermodulation generator of value H j2 o 8V02 shunts a 
capacitor of value Qo and for 8i s an intermodulation generator Ih o 
8\ps2 is in parallel with the linear generator Di 8^2, where both D\ 
and Di are evaluated at o>i + W2. The second-order Volterra kernels 
require the values of the circuit node voltages at coi, u>2, o>3 computed 
in steps (i) and (ii). 

(iv) Solve for the response at toi — 03 and C02 — 03, as in step (Hi). 

(v) Solve for the response at wi + C02 — W3. The second-order 
intermodulation generators used in steps (Hi) and (iv) are replaced by 
third-order intermodulation generators (see Fig. 9b). Each third-order 
generator consists of a pure third-order term plus three second-order 
interaction terms computed from the responses at toi + 0)2, wi — W3, 

W2 _ W3, «l, W2, U3. 

There are two kinds of contributors to the third-order generator 
currents. The first kind expresses the interaction of a second-order 
product with a fundamental signal to produce a third-order product 
[eq. (8)]. The second kind expresses the interaction of three funda- 
mental signals to form a third-order product [eq. (9)]. In terms of 

*v o2 




(a) 



6V o3 



Cox 



:c d0 6L 



Fig. 9 — Equivalent circuits for distortion analysis, (a) Circuit for finding second-order 
voltages, (b) Circuit for finding third-order voltages. 
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the output voltages in volts across load resistor R L , M 3 becomes 

,1 fiVo3(wi + co 2 - co 3 ) I 2 X 10~ 3 R L 

M3 = 2 ° l0g \XV ( \WRV ( l\\RV ( \\ - 15 - 6 - (38) 

1 Vbi(coi) II o V01IW2) || oVoi(ws) I 

The output voltages found in steps (£), {ii), and (u) above are inserted 
into (38) to find the value of M3 . 

IV. DETERMINING SURFACE PARAMETERS 

Several parameters are used in this model to characterize the oxide- 
semiconductor interface. These describe: 

(i) The distribution of surface states across the semiconductor band- 
gap. 

(ii) The capacitance per unit area of the oxide. 

(Hi) The impurity concentration at the surface of the semiconduc- 
tor. 

(iv) The electron capture coefficient of the surface states. 

(u) The fixed charge density in the oxide. 

Although for these parameters the literature provides order-of-mag- 
nitude values that are representative of this kind of interface, some 
values are strongly dependent on the specific processing used to grow 
the oxide. This is particularly true of the surface-state distribution. 
Techniques for measuring surface-state distribution can be found in 
the literature. The surface-state density can be inferred, for example, 
from measurements of the ac conductance of an mos capacitor. 6 

Measurements of capacitance and conductance were available for 
mos capacitors fabricated with the identical oxidation process used for 
the transistors. These data covered wide ranges of frequency (0.1, 1, 5, 
10, and 30 MHz) and bias voltage (from volts into inversion). If 
appropriate values of the parameters were used, the ac and dc surface 
model developed for the intermodulation analysis should be able to 
predict these measurements quite accurately. An optimization program 
has been used to adjust the parameter values in this model to make it 
fit the measured admittance data. The surface-state density has been 
represented by a sixth- degree power series and each of the other 
parameters by a single average value. 

Two small additions had to be made to the equivalent circuit of Fig. 
7 to achieve a good fit of the data. A conductance was added in shunt 
with the oxide capacitance, and a series resistance was inserted at one 
of the terminals. 

The fitting of the data required the use of an adaptive nonlinear 
least-squares algorithm. 11 A set of 12 parameters, seven of which 
characterized the surface-state distribution, were optimized by the 
algorithm. A thirteenth parameter, the electron capture cross section, 

INTERMODULATION DISTORTION IN UNDERSEA CABLE 2039 



was held fixed at 5 X 10~ 16 cm 2 . Figures 10 and 11 show the fitting of 
the conductance and capacitance data for the (111) -oriented capaci- 
tors by the model. Figures 12 and 13 show the fitting for ( 100) -oriented 
capacitors. The surface-state distributions resulting from this data 
fitting appear in Fig. 14; the rest of the parameters are listed in Table 
I. The surface-state distributions of Fig. 14 appear somewhat higher 
and have narrower valleys than is usually found for Si0 2 -Si interfaces. 
It is known that the deposition of SiN layers over the Si0 2 layers can 
modify the surface-state distributions in this fashion. 12 

The values of oxide conductance in Table I are quite consistent with 
recently reported oxide conductance measurements of T. M. Nazar. 13 
Although the frequency range of Nazar's measurements was below 
ours, his measurements show the conductance to be roughly propor- 
tional to frequency over a wide range. From this behavior, he suggests 
that the conductance is due to a hopping process of conduction. 
However, no mention was made of the orientation of his silicon 
substrates. In our model, it would probably have been better to 
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Fig. 10— Conductance of mos capacitor. 
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Fig. 11— Capacitance of mos capacitor. 
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Fig. 12— Conductance of mos capacitor. 
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Fig. 13 — Capacitance of mos capacitor. 

represent the conductance by a simple frequency dependence than to 
assume it to be constant, as we did. In any case, this has little effect on 
the computed modulation coefficient. 

V. DISCUSSION OF RESULTS OBTAINED WITH THE COLLECTOR-BASE 
JUNCTION-CONTROLLED MOS CAPACITOR MODEL 

Using the algorithms described above, values of Ma have been 
computed for two collector-base junctions, each with a different junc- 
tion-controlled mos capacitor. These capacitors are characterized by 
the parameters given in Table I, which are representative of two 
different silicon crystal orientations, (111) and <100), with high and 
low surface-state densities. The collector doping profile of the SG 
transistor is far from flat near the collector-base junction. Hence, to 
find the junction capacitance, we preferred to use an algorithm that 
could calculate junction capacitance for an arbitrary doping profile. 
H.J.J. De Man has described a fast algorithm for such a calculation. 14 
The algorithm used in this model for the collector-base junction 
capacitance is based on De Man's work. 

Rd turned out to be more than two orders of magnitude smaller than 
R a for these simulations and so was ignored. Moreover, in (22) Sxpm 
was assumed to be negligible in comparison with 8\p s . 

Computed and measured values of Mz as a function of collector-base 
bias voltage are plotted in Figs. 15 and 16. Figure 15 represents the 
case of (111) crystal orientation, and Fig. 16 represents (100). The 
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Fig. 14 — Distribution of surface states. 



Table I — Values of surface parameters 



Parameter 


Orientation 


(111) (100) 


Oxide capacitance per unit area, nF cm -2 
Surface impurity concentration, 10 l5 cm -3 
Fixed charge density, 10"/cm 2 
Oxide conductance, micromhos/cm 2 
Series resistance, ohms 


3.83 4.19 

1.08 4.67 

5.85 3.15 

75. 0.32 

17. 1.6 



computed curve for the (111) junction agrees qualitatively with the 
measured data in that it shows the hump that is characteristic of ( 111 ) 
junctions. The appearance of the hump at higher voltage in the 
computed curve than in the measured data is, we believe, due to 
inaccuracies in determining the surface parameters used for the sim- 
ulation and to approximations used in the modeling. 

The position and height of the hump are sensitive to the values of 
fixed charge and average surface-state density used in the simulation, 
as shown in Fig. 17. Even though the mos capacitors used to determine 
the surface parameters all came from a single (111) -oriented silicon 
slice and a single (100) slice, there was considerable variation in the 
measured admittance data for the (111) samples. This, in turn, pro- 
duces a spread in the surface-state parameters derived from fitting the 
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Fig. 15— Intermodulation characteristic. 











<100> ORIENTATION 












28.33, 29.21, AND 29.65 MHz 




x \ 












X 

\ 


V 








- 






V 












\ 


\\ 


X X MEASURED DATA 


- 








\\ 

\\ 

^x\ 

I I I \ 






1 




i 


1 III 



6 8 10 20 

BIAS VOLTAGE IN VOLTS 



60 80 100 



Fig. 16— Intermodulation characteristic. 

admittance data. And this spread undoubtedly contributes to the 
difference between the computed and measured M$. 

As discussed earlier, the dominant surface effect on intermodulation 
would be expected to arise from the region under the base connecting 
pad. And so we have restricted our surface modeling to this portion of 
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Fig. 17 — Effects of changing surface parameters on intermodulation characteristic. 

the intersection of the collector-base junction with the surface. How- 
ever, there are probably also contributions from the extensive region 
not under the base pad where the collector-base junction intersects 
the surface. This region would be much more difficult to model, and 
we have not attempted to do so. This omission is probably responsible 
for some of the difference in the computed and measured Ma, also. 

The computed Ma for the (100) orientation agrees well with mea- 
sured data. Note the almost complete absence of a hump in these 
curves. This is because of the low surface-state and fixed charge 
densities for this orientation. 

If, in computing Ma, for the ( 111 > case, the surface-state nonlinearity 
is omitted, the computed Ma comes out to be that shown by the dashed 
line in Fig. 15. This shows clearly that the surface states are the 
primary source of the difference in intermodulation produced by the 
(100) and (111) -oriented transistors. 

Contour maps of Ma in the fi, fa frequency plane for a fixed value of 
the product frequency (/i + fi — fa) reveal those combinations of 
frequency that give rise to high intermodulation. These maps have 
been computed using the collector-base junction-controlled mos ca- 
pacitor model. Figures 18 and 19 show Ma maps for the (111) junctions 
at bias voltages corresponding to the peaks of the characteristic Ma- 
bias humps (see Fig. 15). These bias points (17 V for the simulation 
and 5 V for the measurements) have been chosen because surface 
effects dominate the distortion at these points. The predicted levels of 
Ma and the general shape of the contours agree reasonably with 
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Fig. 18 — Computed contour map of modulation coefficient. 

measurements. However, the measurements show a stronger frequency 
dependence. 

Contour maps of Ms for (100) junctions are shown in Fig. 20 for the 
model and Fig. 21 for measured data. In the absence of any hump in 
the Ms-bias curves, a bias voltage of 5 V was chosen for the comparison. 
Again the predicted map agrees in magnitude and general shape with 
the measured map. The M 3 is seen to be less frequency-dependent 
than for the (111) case, except at frequencies very close to the product 
frequency. 

It should be pointed out that the nonlinear model only includes 
terms up to third order. This approximation is valid only if the power 
in the fundamental driving tones is sufficiently low. The measured 
data shown in Figs. 15, 16, 19, and 21 have been measured with power 
in each fundamental P set at -5 dBm. Lower power was not used 
because the resulting distortion power is near the test set noise limit 
of -140 dBm when the bias voltage is near 5 V. The power dependence 
of Ms is shown in Fig. 22 for both <111> and <100) crystal orientations 
for a bias voltage of 5 V. 
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Fig. 19— Measured M 3 frequency map for (111) collector-base junction. 

The noticeable rise in M 3 as P is reduced indicates that the distor- 
tion mechanism is not accurately described by only second- and third- 
order terms when P Q is greater than -15 dBm. This is because a pure 
third-order nonlinearity does not yield a power-dependent M 3 . For this 
reason, the model of the collector-base junction does not predict M 3 
variation with power. The lack of higher-order terms in the model may 
also explain the low predicted sensitivity of M 3 to the choice of 
fundamental frequencies. 

The use of Volterra series and the perturbation method to describe 
a complex nonlinear system becomes rather unwieldly when terms 
higher than third-order are included. A simpler analytical method is 
needed to simulate the higher-order nonlinear effects in a collector- 
base junction-controlled mos capacitor. Simulations using the present 
model suggest that the dominant source of surface-state distortion is 
due to the exponential dependence of surface electron density n, on 
surface potential \p s as given in eq. (20). The G 2 and G 3 Volterra kernels 
are negligible contributors to the D h D 2 , and D 3 composite Volterra 
kernels. If a higher-order nonlinear model were to be developed for 
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Fig. 20— Computed contour map of modulation coefficient. 

surface-state intermodulation, perhaps only the n s nonlinearity would 
need to be included. 

These computations have shown how important the contribution of 
surface states can be to the nonlinear performance of collector-base 
junctions. If intermodulation distortion is critical in the application of 
these junctions, it is important to control the distribution of surface 
states at the Si-Si0 2 interface. This is not only a matter of the choice 
of crystal orientation, but a matter of the oxidation process used and 
also of post-oxidation processing. 

VI. MODELING INTERMODULATION DISTORTION IN THE COMPLETE SG 
TRANSISTOR 

Transistor models used in circuit analysis programs such as the 
Gummel-Poon model 15 and Extended Ebers-Moll model 16 are adequate 
for low-current modeling where temperature effects, base pushout, and 
lateral spreading are not important. Low-current operation means 
that, for a fixed collector-emitter voltage, the collector current is below 
the currents for which fa or /? peak. High-frequency circuits require 
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Fig. 21 — Measured Ma frequency map for (100) collector-base junction. 

the highest f T possible for the critical transistors, so the transistors 
must be biased with high currents. 

The high current model presented in this paper extends the Gum- 
mel-Poon model by defining two base charge terms. One is the actual 
base charge and the other is a pseudo-base charge that allows Gum- 
mel's equation to describe the dominant component of collector cur- 
rent even when two-dimensional effects are important. In the limiting 
case of a one-dimensional structure, both charges are identical. 

A number of high-current base-charge models have been pro- 
posed, 17 " 20 but they are based on device physics and need device 
geometry and material parameters as inputs. These models handle 
one-dimensional high current effects accurately, but only approximate 
two-dimensional effects. The base charge and the pseudo-base charge 
discussed in this paper are not derived from device physics but use 
empirical formulas as first suggested by Choma. 21 The equations for 
base charge given by Choma are found to model the quasi-base charge 
adequately, but a new set of equations is used to model the total base 
charge. 
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Model parameter extraction is done with a computer program that 
uses measured currents and voltages of the test transistor as an input. 
A two-pass operation is needed to get a complete model. The first pass 
computes low current model parameters and plots the base charge 
multiplication factor. The user must define a functional form that can 
approximate the base charge multiplication factor as a function of 
collector current and collector-base voltage. The functional form may 
differ between transistor types due to different doping and structure. 
The user then computes initial guesses at the model parameters. The 
second pass uses an optimization program to refine the approximate 
model parameters. The use of a model parameter extraction program 
insures unique model parameters because the parameters are deter- 
mined in precisely the same order with the same assumptions for every 
extraction run. 

6. 1 Transistor model description 

The cross section of a typical bipolar device is shown in Fig. 23. Two 
components of dominant collector current I ecu Ices and base charge 
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Fig. 23 — Transistor cross section. 



Qbi, Qbs associated with the intrinsic region directly under the emitter 
and the sidewall region are defined. The intrinsic region has current 
flow in only one dimension, and Gummel's equation applies where 



Icci = 



IsQBo[exp(VWVr) - exp(V BC /V r )] 
Qbi 



(39) 



Of course, this assumes that lateral base current flow has a negligible 
effect. Current flow in the sidewall region is two-dimensional, and no 
rigorous analog of Gummel's equation exists to describe Ices- However, 
if one sacrifices rigor, (39) can be modified to approximate Ice 



t t ^r IsQBo[ex p(V B E/VT) - «p(Vbc/Vt)] 

■ice — icci + Ices = — , 

Qbq 



(40) 



where Qbq is defined as the quasi-base charge which equals Qbi for a 
purely intrinsic device. At low collector currents, Qbq « Qbi. However, 
at high currents both base pushout and lateral spreading occur, and 
Qbq must change; so (40) models Icci + Ices. Qbq can be expressed as 
shown below. 



Qbq — Qbo 



1 + 



Qi 



CjeoVbo 



+ n Q v + 7T-B x (Vbe, V bc )Icc 



(41) 



Qbo = approximate zero bias total majority base charge 
Qe = base-emitter junction majority charge 
Qc = base-collector junction majority charge 

Vao = forward Early voltage 

Vflo = reverse Early voltage 

Tpo = extrapolated low-current transit time for Vcb = 
Bx(Vbe, Vbc) = base pushout and lateral spreading function. 

For an intrinsic device, CjeoVbo = CjcoVao = Qbo. However, this is not 
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true for a real two-dimensional device. The emitter junction charge 
Qe and the associated capacitance C E are given below. 



Qi 



rV BE 

Jo 



C E (V) dV 



(42) 



C E = 



C Ea »x Vbe>Pe(1-B e ) 

-m s (43) 



yJEO 



Vbe 
Pe 



Vbe < Pe(X - B E ) 



Be 



\CEmax/ 



(44) 



The zero-bias capacitance is Cjeo, whereas the peak value of the 
depletion capacitance in forward bias is G&m»- Be defines the voltage 
at which C E becomes constant at C E max, and P E is the built-in voltage. 
The collector charge and capacitance are defined with similar equa- 
tions and model parameters Cjco, Pe, M c . An empirical function for 
B x is described by Choma 21 and is given below: 

B x (Vbe, Vbc) = Bxo + (1 - Bxo)e- f (45) 

B X o=l + ae v ™ /b (46) 



f=c(e 



VBE/d _ 



1), 



(47) 



where a, b, c, and d are model parameters. 

This function is adequate for test devices studied in this paper, but 
other functions may be necessary for different device structures. 

Equation (41) defines Qbq in terms of Ice', to define Qbq in terms of 
only Vbe and Vbc, substitute (40) into (41) and solve for Qbq to get 

Qbq a l/ H Qe | Qc 
Qbq 2 V CjeoVbo CjcoVao 



1 
+ 2 



1 + 



Qi 



Qc 



Cjbo Vbo Cjco Vao 



+ — — axis 



Q 



'BO 



€ Vbe/V t _ qVbc/Vt ) 



-i 1/2 



(48) 



The base-pushout and lateral spreading function Bx as defined by (45), 
(46), and (47) is only valid for the active or quasi-saturation regions of 
operation. A more complete description of Bx is needed to model 
inverse operation. 
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The description of Ice given above is combined with the other 
elements of the Gummel-Poon model 15 to describe base current and 
charge storage as shown in Fig. 24. The base currents injected into the 
collector Ibc and emitter Ibe and the base resistance Rb are given by 



Ibe = he v °* /VT + I* V **' N * V * 

I BC = I 3 e v *c/VT + l4e v BC /N c v T 

Rb = Rbx + Rbi/(Qbq/Qbo)- 



(49) 
(50) 
(51) 



The base resistance due to the external base is modeled with a constant 
resistance Rbx and a conductivity-modulated resistance Rbi. 

Self-heating of a device is an important consideration for high- 
current modeling, so an auxiliary model is included. A current gener- 
ator equal to the power being dissipated in the device drives the 
thermal resistance Rt and capacitance Ct- The voltage across Rt 
equals the instantaneous temperature of the device. For low-frequency 
circuit simulations, the temperature of the device must be included as 
a dynamic variable in all the model equations. However, for high 
frequency simulations, only the dc operating point must be computed 
with temperature as a variable, since the time constant of the RtCt 
network is long. 





Ct 



Fig. 24— Transistor model topology. 
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Charge storage is modeled by capacitors Cbe and Cbc where 

r _ *Qb p _ dQ B 

Cbe = ~^—, Cbc - ~rz- (5^) 

a Vbe o Vbc 

and Qb is the actual charge in the base not to be confused with Qbq. 
This definition of Cbe and C B c assumes the quasi-static approxima- 
tion 22 for charge distributions in the base. The base charge distribution 
at time fo is given by the dc charge distribution evaluated for Vbe = 
Vbe (to) and Vbc = Vbc (to)- The excess phase in the ft versus frequency 
measurement is caused by delay in the base and in the collector space- 
charge regions. The quasi-static approximation does not model the 
diffusion of carriers in the base or drift of carriers in the space-charge 
region; hence, excess phase is not modeled. However, at frequencies 
that are low compared to the fa of the device, the quasi-static model 
is useful. 

The actual base charge consists of a constant charge due to doping 
Qbo, depletion charge Qe, Qc, and diffusion charge, and is defined in 
a similar fashion as Qbq to be 



Qb = Qbo + Qe + Qc + TpoBy(Icc, Vbc)Icc- 



(53) 



Note that Qbo is not really necessary in (53) since only the derivatives 
of Q B are used. The function By (Vbe, Vbc) accounts for the effects of 
base pushout and lateral spreading. It is necessary to use a different 
functional form for By than used for B x . The fa vs Ic data cannot be 
matched closely at both low and high V c if By has the same form as 
B x . However, with By given by the function in (54), the fa data could 
be matched over the entire range of Vcb as shown in Fig. 25. The dc 
/? vs Ic curve corresponding to Fig. 25 is shown in Fig. 26. 

By = Alec + Byo (54a) 

n - 1 + me- Ic/ni (54b) 

Byo - B ie VBc/B > + B 3 (54c) 

A = AxV%c + A 2 V B c + A 3 (54d) 

no = Ci + C 2 Vbc (54e) 

m = C 3 + C 4 Vbc (54f) 

with model parameters A\, A 2 , A 3 , Bi, B 2 , B 3 , d, C 2 , C 3 , C 4 . Equation 
(54) represents a function of considerable complexity. However, each 
subfunction (54b) to (54f) is determined from plots of A, n, and Byo vs 
Vbc obtained from measured data. The exact forms of these equations 
may differ for devices with completely different structures. 
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VII. COMPUTATIONS USING THE TRANSISTOR MODEL 

The SG transistor has been the product of years of effort to develop 
a highly linear transistor for submarine cable use. 2,3,24 As a result, 
extreme flatness of the Ibf. versus Vbe and /? versus Vnc and I c 
characteristics has been achieved."' 24 

We have tried to model the SG transistor carefully, so that the 
model would reproduce measured distortion accurately when anoma- 
lous distortion is absent. However, in this attempt, it has become 
apparent that the distortion predicted by the transistor model per se 
without surface effects has been excessive. 

The primary source of the excess transistor distortion from the 
model is its inability to adequately simulate the extreme linearity of 
the SG transistor. To accomplish this involves fitting not only fa and 
/? accurately over a wide range of voltage and current bias, but the 
second and third derivatives of these quantities as well. Poon 24 has 
shown that, in the high frequency, limit M 3 is strongly related to the 
second derivative of fr. At lower frequencies, ft becomes important. 
Unfortunately, the modeling equations are not capable of matching /3 
and fr over a wide bias range. To improve the model, the functions 
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Fig. 26 — Direct-current beta-current characteristics. 

describing base charge Bx (47) and B Y (54) must be increased in 
complexity so as to fit the higher derivatives of charge. 

Replacing the collector-base diode with the full transistor model 
produced no fundamentally new problems to the intermodulation 
computation itself. The principal effect was simply to increase the 
complexity of the computations as a result of the increase in equivalent 
circuit complexity. 

Intermodulation analysis requires derivatives up to third order for 
all nonlinear elements. These are easily obtained explicitly except for 
the normalized base charge Qb, which is a function of Vbe and Vbc 
Derivatives of Qb have been computed numerically to avoid derivation 
mistakes and to allow quick changes in the equations for Qb without 
rederiving the derivatives. 

For an operating bias of Vcb = 10 V, the internal junction temper- 
ature is elevated about 75 °C. The inclusion of junction heating in the 
transistor model lowers the simulated Ma by some 4 dB at the higher 
voltages. This occurs because of the slight reduction of exponential 
modeling factors of the form e gV/kT . 
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Computed and measured plots of M% versus collector-base bias are 
shown in Fig. 27. The measured data show the characteristic hump for 
the < 111) -oriented transistor. These plots are very similar to those for 
just the collector-base junction (Figs. 15 and 16), the principal differ- 
ence being that the transistor amplifies the distortion by about 10 dB. 

For the simulated Mz curve, the fi transistor nonlinearities were 
artificially reduced in accordance with the comments above. Without 
this reduction, the M 3 curve from the transistor alone would be about 
14 dB higher at the high-voltage end and would override any contri- 
bution from surface effects. The effect of the surface-state contribution 
to Ma is seen to be a hump at approximately the same level as the 
measured hump. 

Although the transistor model does not allow exact prediction of 
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Fig. 27— Transistor intermodulation-voltage characteristics. 
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distortion, simulations were made to show variations of Mz with 
collector current. As shown in Fig. 28, these follow the same trend as 
measured data. The general shape of the curves is a result of transistor 
nonlinearities, not surface effects. Base current is composed of two 
terms. 



T IcdQ B 

lB = J+-dT- 



(55) 



At high currents where /? is constant, Qb is dominated by depletion 
charge, which is nonlinearly related to Ic. As Ic increases, Qb becomes 
dominated by diffusion charge, and Qb is proportional to Ic/fr. The 
first minimum in Mz is reached at the current for which dfr/dlc = 0, 
that is, at the peak of the fr — Ic curve. The next local minimum in M3 
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Fig. 28 — Transistor intermodulation-current characteristics. 
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occurs where f^fr/dl* = 0, that is, at the inflection point in the fa — Ic 
curve. This simplified discussion is valid at high frequency where 
dQa/dt dominates in (55). For lower frequencies, the effect of nonlin- 
earities in /? becomes important. 
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APPENDIX A 

Integrating Volterra Kernels over the Surface-State Distribution 

As can be seen from equation (19) and Appendix C, to find the total 
current flowing into the surface states involves a sum of integrals 
containing the factor 

*»-» (56, 



1 + jufo/CnTlso ' 




where f = [1 + g s exp(u — Mf)]" 1 . 




Now 




f(1 f) - 8sexp(u-u F ) 


df 


/oU Jo) r . 

[1 + g s exp(u - u F 


l 2 du 



(57) 

is sharply peaked about u = uf. This is readily seen by plotting fo(u) 
and noticing that it slopes steeply only near u = uf. Since dfo/du is 
peaked, so are fo(df /du) and fl (dfo/du). If w/c„/iso <c 1, that is, at low 
frequency and high semiconductor doping, (56) reduces approximately 
to (57). Therefore, provided that N a is essentially constant through 
this peak, the integrals in (19) can be readily evaluated in the manner 
of Nicollian and Goetzberger. 7 

However, as w/c/iso becomes comparable with or greater than unity, 
i.e., at high frequency or low doping, the peak broadens and shifts 
away from uf toward the conduction band. The location of the peak 
can be found by setting the derivative of (56) with respect to u equal 
to zero and solving for u. The integrals could then be evaluated as 
above with N s taken out of the integral and set equal to the value at 
the peak location. But broadening the peak makes this approach rather 
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unpalatable, and we have preferred to evaluate the integrals numeri- 
cally. 

APPENDIX B 

Taking into Account Surface Potential Fluctuations in Calculating Surface- 
State Volterra Kernels 

Following Nicollian and Goetzberger, 7 assume that the dc surface 
potential is uniform over a small region of characteristic area a. The 
probability that the dc surface potential on a particular one of these 
regions lies between u 8 and u 8 + du 8 is P(u s ) du 8 . Then, to find the 
first-order Volterra kernel A(wi), we multiply (23) by P(u a ) du 8 and 
integrate over all possible surface potentials. 

f°° f°° a f" c 

Di(ui)P(u s ) du 8 = I ywi— I N s (u, u s )Hi(o}i) du 

J—aa J— oo Jlly 

• nsoP(u 8 ) du„. (58) 

We recognize that /o is not a function of u 8 and assume P(u 8 ) to be 
sharply peaked at the average surface potential u„. In particular, we 
assume P(u 8 ) to be Gaussian. 

-{u 8 - u 8 ) 2 ' 



P(u 8 ) = (2™ 2 )~ 1/2 exp 

where 

A 
o 8 = 



2oi 



(59) 



qQ 



-1 1/2 



a(D d ) 



(C ox + Gi ) V T 

and Q is the mean oxide charge density in C/cm 2 . 

Then inverting the order of integration, we assume that N 8 {u, u 8 ) 
varies slowly with u s and can be taken outside the integral with respect 
to Us. We then have 



I 



Di(ui)P(u 8 ) du s 



•J Uy ** 00 

where u 8 is the average value of the surface potential and 

JCOifo 

R(u s ) = exp [— (u s — Uf)\. 

CnTli 

Now expand [1 + R(u 8 )Y l in a Taylor's series about u 8 , 

1 F"(u) 

F ^ U ") = , . n/ =F(u 8 )+F'(u 8 )(u 8 -u 8 ) + - (u 8 - u 8 ) 2 + .... 

1 -I- R{u 8 ) 2 
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Then 



. , * , du s « F(u s ) P(u 8 ) du 8 

+ F'(u.) I (n.-iZ.)P( M .)rfifc 
+ — g - " (" 8 ~ u a ) 2 P(u s ) dus 



*F(Us)+-^F"(a 8 ). 



The second integral on the right vanishes. Hence, 



(61) 



£ 



A(coi)P(ii.) du s = 



JuiQ 



u c 



Nsfoil - fo)Pi(w h Us) du, (62) 



Uy 



where A(wi> w 8 ) is the quantity on the right side of (61). 

Treating the second- and third-order Volterra kernels similarly, we 
obtain 

J A(6Ji, (o 2 )P(w 8 ) rfu 8 s - ^t^ j C [^( Wl + <o 2 , £)/„ 

J-ee T Ju v 

- A(wi, (0 2 , M.)/g]iV.(ll)(l - /o) cfa, (63) 
where, letting (Oi + a>2 — (03 = (o p , 

Fa(<Ol, (02, M s ) = F((Oi, COi + (02, M s ) + P((02, (Oi + (0 2 , M fl ) 



+ ~T [-^"("1. "1 + "2, COp, W 8 ) 
+ F"((0 2 , Wl + (02, (Op, M s )] 



£ 



A(<oi, co 2 , (o 3 )P(u 8 ) cfa s s -^~2 



-^l(Wp, Ms) 



/o 



— i^2(<0i, (02 + (03, W 8 )/o 

P 3 (<Oi, (0 2 , (0 3 , Um) . 
H ~ /< 



■ 



■(l-/o)iV.( M )dtt, 



(64) 
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where the bar is used as in (25) and 

1*3(W1, t*>2, "s) = F( U h W l + U 2, 0)p, U«) + F(b)2, 0)i + 0)2, 0)p, U a ) 

+ F(o)i, 0)1 + 0) 3 , 0)p, M s ) + F(o) 3 , Wi + C03, 0)p, Us) 
+ F(032, 0)2 + 0)3, 0)p, Us) + F(0)3, 0)2 + 0) 3 , 0)p, U,) 

+ y [-F"(o)i, Wl + 0)2, «p, Us) 

+ F"(0)2, 0)1 + 0)2, 0)p, Us) 

+ F"(0)i, 0)1 + 0)3, 0)p, Us) + F"(0)3, 0)1 + 0)3, 0)p, U 8 ) 

+ F"(0)i, 0)2 + 0)3, 0)p, Us) + ^"(0)3, 0)2 + 0) 3 , 0)p, U s )]. 

The definitions of F and F" are as follows: 

1 



F(<Ol, 0)2, • • • , w n ) = II 



ti l + RM 



F"(U)1, 0)2, • • ', 0)„) = F(0)l, 0)2, • • •, 0)„) 



I R(ui)F(u t ) - 2 R(ui)F 2 ((*) 

i=l / i-l 



APPENDIX C 

Derivation of Surface-State Vo/terra Kerne/s 

A Volterra series for / in terms of n s , eq. (12), describing the time 
dependence of the Fermi level is rewritten as 



— = c„(l - f)n a - e„f. 
at 

First break (56) into a dc part and an ac part. Let 

f = fo + 8f, n s = fiso + Sn s 

4-Jf= c n (l - fi> - Sf)(nso + 8n 3 ) - e n (fo + «/). 
at 

Direct current equation: 

= c„(l - fo)ngo - enfo. 



Therefore 



or 



e„ = • 



C n (l — /o)/lgQ 

/o 



CnflsO 



(65) 



= Cnfiso + e n . 
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Alternating current equation: 



— S/= c„(l - f )8n s - (c n n s -o + e n )Sf- CnSfSn* 
at 



or 



£-8f= c„(l - f )8n. - ^ 8f- CnSfSn,. (66) 

a£ /o 

It is desirable to express Sf as a Volterra series 

8f = Hi(a3 a ) o 8n s + 7/2(0)0, to/,) o Sni; + Hz(u a , to/,, to c ) o 5/i'j, (67) 

where the o operator indicates that the magnitude and phase of each 
term in 8n" is to be changed by the magnitude and phase of iJ„(to„, 
• • •, to,,). To find the Volterra kernels A„{u\, • • -, to,,) substitute (58) 
into the ac eq. (57) using a theorem for higher order Fourier transforms 
of derivatives (76) to get: 

ju a H\(Ua) O Sn s + j(u) a + WfrJfl^Ua, «*) O 5/1* 

+ y(to a + tO/, + U c )H:)(u a , tOft, to r ) o 8ni 

= c„(l - /o)5n s — ^A [7/i(to a ) o 6n.s + H 2 (u a , to/,) o <5n; 

+ HziUa, 10/,, C0 f ) O 5/lf] 

— c„{[//i(co a ) o 8n s ]8n s + [//2(to a , to/,) o 8n%]8n a 
+ [H:i(u) a , to/,, Uc) o 8n.s]8n s ). 

Consider 8n s = a,\ cos(toi£ + <£i) + a? cos(to 2 ^ + 4)2) + fla cos(to,^ + <>.t) as 
the forcing function and ignore the [H:i(u a , co&, co c ) o 8nl]8n b term since 
it contributes fourth-order terms and small second-order terms. Also 
ignore the small contribution to the first-order terms due to [H>{u tl , 
to/,) o 8ni]8n s . Equate terms of the same order to get three equations: 

Ju„Hi(b)a) o 8n s = c„(l - f n ) o 8n„ — ^— HAu a ) o 8n* (68) 

/<> 

y(to„ + u)b)H-2(u a , to/,) o 5n: = — H2i.Ua, to/,) o 8m 

/n 

— c,[//i(tOa) o 6n.s]6n.s (69) 

_/(t0„ + CO/, + (t) c )H:i(ui a , COfe, tO,) O 5/1* = -— Ha(<*)a, tOfc, tO, ) O 5aI"* 

/o 
- c n [Ho(to„, to/,) o 5n;]5n s . (70) 
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Solving (68) for Hi gives 



HiitOa) =- 



Cnd ~ fo) 



J03 a + (CnJlgo/fo) 

Solving (69) for H 2 using (77) gives 

j(u3 a + o3 b )H 2 (o3 a , Mb) o8n 2 s = — ^— H 2 (u3 a , ojb) o 8nl 

to 



--^[#i(co a )+iJ,(w 6 )]oSn 2 . 



So 



Cntif^Wa) + Hi(m)] 
M2K.Ua, 03 b ) = p 



j(o3 a + 03 b ) + 



CnriM 



fo 



Solving (70) for H 3 using (78) gives 

j(o3 a + 03b + 0} c )Ha(0)a, 03b, 03c) O 5/lg = — Hz{o3 a , 03b, 03 c ) o 8n a 

to 



So 



H 3 (03 a , 03 b , 03 c ) = ~ 



— ■=■ [H 2 (03 a , 03b) + H 2 (03 a , 03 c ) 
O 

+ H 2 (ojb, o3 c )] o Sri*. 

C n [H 2 (03 a , 03b) + H 2 (03 a , 03 c ) + H 2 (o3 b , 03 c )] 



j(03 a + 03b + 03c) + 



CnflsO 



fo 



APPENDIX D 

Using Volterra Series 

Given a nonlinear system with an input u(t) and output y(t), one 
may describe the system using a functional series of the form 2 

y(t)=yAt)+y 2 (t) + --yn(t) t (71) 

where 

yn(t) = I • • • I h n (t - Ti, • • •, t - T„)«(ti), • • •, U{r n ) dTa, '•', d.T n . 

The yi(t) term is the ordinary convolution integral that is used for 
linear system analysis and hi{t) is the impulse response. The kernel 
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h n (t) is the general impulse response. If the input u(t) is changed by a 
gain factor e, then the new output is 



y>n(t) = | •■• I h n (t- Ti, ■■■,t-T n )eu(Tl),---, 

€u(r n ) dr u • • •, dr„ = e n y n (t). 

The system can be thought of as a parallel combination of first, second, 
third, • • • nth order subsystems. 

The following theorems deal with Volterra series in the frequency 
domain (7) to (9) using the circle operator defined in (17). 
Let 

U = A(cOi) O COS U\t + • • • + A(ua) O COS Oiat, 

then 

Hi(co a ) ou = Hi(ui)A(ui) o cos uit + H\(u 2 )A(u 2 ) o cos w 2 t 

+ Hi(w3)A(ua) o cos uat 

= Hi(u a )A(u a ) ocos u a t a = 1, 2, 3 (72) 

H 2 (0J a , Ub) o u 2 

= dc terms 

+ £2(0)1, ui)A 2 (o3i) [Vfc] o cos 2toit 

+ H 2 (oj 2 , u 2 )A 2 (u 2 ) \}/z\ o cos 2u 2 t 

+ H 2 (o)3, u>3)A 2 (ua) [ ] /£] o cos 2u>Wat 

+ H 2 (coi, — (o 2 )A(ui)A(— u 2 ) o cos(wi — u 2 )t 

+ H 2 (it3\, oj 2 )A{loi)A(co 2 ) o cos(wi + u 2 )t 

= H 2 (co a , Ub)A(u a )A(u b ) [%] o cos(o> a + u b )t, (73) 

where [V2] is used for dc and second-order harmonics only. 

Ha(u a , Ub, 00c) O U 3 

= lst-order terms 

+ £3(0)1, coi, u\)A 3 (u)\) [%] o cos 3toi£ 

+ Ha(u 2 , W2, C02)i4 3 (cO>) [WJ O COS 3t02< 

+ £3(0)3, (1)3, W3)A 3 ((0a) [%] o cos 3o>3£ 

+ £3(0)1, 0)1, — W2)A 2 (cOi)A(— CJ2) [%] o cos(2wi — C02)^ 

+ £3(0)1, 0)1, (02)^4 2 (wi)i4(u)2) [%] o cos(2coi + (02)^ 
: (74) 

+ £3(0)1, W2, aj3)A(t0i)A(aJ2)A(aj 3 ) [%] O COS(o)i + a>2 + WsK 

-I- £3(0)1, —0)2, U)3)A(Wi)A(— 0)2)^(0)3) [%] o cos(coi — CO2 + 0)3) f 

=s H 3 (co a , cob, tu<.)A(<Oa)^4(t06)A(<o c ) [Fraction] o cos(to a 4- ctb + u c )t. 
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This last formula is useful for all the third-order terms but not the 
first-order terms. 

[HAlOa) O Uf = %[HAWa)H 2 (Ub, Wr) + Hl(Ub)H 2 (cO a , U c ) 

+ HyMH 2 (o)a, w ft )] o u 3 



= H\MHz(ub, u>,) o u 3 , (75) 

where the bar is used as in (25). 
Differentiating Volterra series, 

—yM < — > (ywj + • • • +/ujr)Hjr(fi>i, ■ ••, wa) (76) 

dt 

[Hi(u„) o«]x« = ai| ifi(wi)| cos(toi< + 0i + <f> u i) 

+ a, | ifi(w 2 )| cos(w-2< + 02 + 0«2) 

+ a:»|Hi(<0.i)| COS(W3< + 03 + <M] x "• 

For harmonic terms, i.e., 
or 

03 a = 60] COft = — <0 a — > dc 
[Hi{03 a ) OU]X U= Vz A(u a )A(0)b)HiM o cos(u„ + U b )t. 

For sum and difference terms, 

[i/i(co a ) o u]x u 

= v h A(o)a)A(u b )[Hi(u a ) + Hibtb)] o cos(w a + Ub)t. 

In general, 

[HiM ou]x u = -AMA(u>b)[HAua) 



+ Hdu b )] 



O COS(<0 a + Ub)t 



= \[HxM+HAub)]ou 2 (77) 

[H 2 {u a , Ub) o u 2 ]X u = % [H 2 (cj a , u b ) + H 2 (io a , u c ) 

+ H 2 (cJ b , We)] o u 3 . (78) 

APPENDIX E 

List of Symbols 

B.v, By Base pushout and lateral spreading functions. 

C,i Surface depletion capacitance. 
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Ci Surface inversion capacitance. 

c„ Electron capture probability. 

Cox Oxide capacitance. 

C a Surface-state capacitance. 

Ct Thermal capacitance. 

d Depletion layer width. 

D n Complete Fourier-transformed surface-state Volterra ker- 
nel of order n. 

E Energy level. 

e„ Electron emission probability. 

fo, f(t) Direct-current and time-dependent fermi distribution 
functions 

G n nth order component of Fourier-transformed surface-state 
Volterra kernel. 

g a Ground state degeneracy. 

h n Time-dependent Volterra kernel of order n. 

H n Fourier-transformed Volterra kernel of order n. 

Ibc Base current injected into the collector. 

Ibe Base current injected into the emitter. 

Icci, Ices Intrinsic and sidewall components of dominant collector 
current. 

i a (t) Time-dependent surface-state current. 

k Boltzmann's constant. 

Ld Intrinsic Debye length. 

Ms Third-order modulation coefficient. 

m Intrinsic electron density. 

N s Surface-state density. 

fiso, n s (t) Direct-current and time-dependent electron surface con- 
centrations. 

Po Output power. 

q Electronic charge. 

Qb Actual charge in base. 

Qbi, Qbs Intrinsic and sidewall components of base charge. 

Qbo Total majority base charge at zero bias. 

Qc Base-collector majority charge. 

Qe Base-emitter majority charge. 

Q g Gate charge. 

Qi Surface inversion charge. 

Q s Semiconductor charge. 

Q ss Surface-state charge. 

Qo Fixed oxide charge. 

Q Mean oxide charge density in C/cm 2 . 

Rb Base resistance. 

Rd Surface depletion layer resistance. 

Ri Surface inversion layer resistance. 
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R nD Surface depletion layer resistance for electron flow. 

Rns Surface-state resistance for electron flow. 

R pB Resistance to hole flow in quasi-neutral region adjacent 

to depletion layer. 

R pD Surface depletion layer resistance for hole flow. 

R ps Surface-state resistance for hole flow. 

R T Thermal resistance. 

t, to Time. 

T Temperature. 

Tfx, Low-current transit time for Vcb = 0. 

u Fn , ufp Fermi potentials for electrons and holes in units of the 

thermal voltage. 

u 8 Surface potential in units of the thermal voltage. 

Vao Forward Early voltage. 

Vbc Base-collector voltage. 

Vbe Base-emitter voltage. 

Vbo Reverse Early voltage. 

V T Thermal voltage (kT/q). 

Y a Surface-state admittance. 

8 Alternating current component of. 

(i n Electron mobility, 

co Angular frequency. 

$Fn Electron quasi-fermi potential in electron volts. 

4>'ms Metal-semiconductor work function. 

\p B Surface potential in volts. 
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